Quantized antiferromagnetic spin waves in the molecular Heisenberg ring CsFeg 
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We report on inelastic neutron scattering (INS) measurements on the molecular spin ring CsFeg, 
in which eight spin-5/2 Fe(III) ions are coupled by nearest-neighbor antiferromagnetic Heisenberg 
interaction. We have recorded INS data on a non-deuterated powder sample up to high energies at 
the time-of-flight spectrometers FOCUS at PSI and MARI at ISIS, which clearly show the excitation 
of spin waves in the ring. Due to the small number of spin sites, the spin-wave dispersion relation 
is not continuous but quantized. Furthermore, the system exhibits a gap between the ground state 
and the first excited state. We have modeled our data using exact diagonalization of a Heisenberg- 
exchange Hamiltonian together with a small single-ion anisotropy term. Due to the molecule's 
symmetry, only two parameters J and D are needed to obtain excellent agreement with the data. 
The results can be well described within the framework of the rotational-band model as well as 
antiferromagnetic spin-wave theories. 

PACS numbers: 75.50.Xx, 75.10.Jm, 78.70.Nx 



I. INTRODUCTION 

In recent years, large progress has been made in the un- 
derstanding and control of mesoscopic magnetic systems, 
which is driven by the interest in fundamental research 
as well as potential applications in molecular devicesii"— 
Dimensionality plays a key role here: magnetic confine- 
ment in one or more dimensions leads to strong changes in 
the density of states, i.e., quantization effects occur j^"— 
which become more pronounced with lower dimensional- 
ity. Theoretical work on spin waves in low-dimensional 
antiferromagnetic (AFM) systems reaches back to the 
early 1950s. While low-temperature thermodynamical 
properties in two and three dimensions are well described 
by a conventional spin-wave approach, the case of one 
dimension (ID) is more difficult due to divergencies in 
the sublattice magnetizationJ^ii^ These divergencies are 
deeply rooted in the fact that, by construction, spin- 
wave theory (SWT) starts out from the assumption of 
a long-range ordered ground state, and their appear- 
ance provides a hint to the nature of the ground state, 
i.e., whether it is disordered (divergency) or long-range 
ordered (no divergency). However, despite these sub- 
tleties it has been found that also for disordered sys- 
tems SWTs can produce results of high accuracy, e.g., 
for the ground-state energy.Ji Furthermore, the divergen- 
cies could subsequently be removed by introducing con- 
straints on the sublattice magnetization in the so-called 
modified SWTsi^-^ and finite-size SWTs^"^ This es- 



tablishes a significant improvement, but nevertheless the 
applicability of SWTs to disordered systems remains a 
topic of current research. Formerly, the main focus of 
the theoretical and experimental activities has been on 
systems of infinite extensionj^iSS, and it has only been in 
recent years that research on finite, or zero-dimensional 
(OD), systems has been intensified. In this context, the 
topic of quantized spin waves in mesoscopic spin clusters 
has been investigated 4*2^""— 

A good physical realization of such mesoscopic spin 
clusters are molecular nanomagnets.~ They have been 
shown to be a vast experimental playground enabling the 
study of precisely defined clusters with chemically engi- 
neered exchange couplings and magnetic anisotropy, or 
topology. So far, a variety of mesoscopic and quantum 
phenomena such as tunneling of the magnetization and 
stepwise change in the magnetization have been demon- 
strated (for reviews see Refs. IssI - IstI ). A particularly in- 
triguing type of spin clusters are AFM Heisenberg rings, 
in which a number of magnetic metal ions are arranged in 
the shape of a ring with nearest-neighbor AFM Heisen- 
berg couplings. The Hamiltonian describing their mag- 
netism is given in its simplest form by 

H = - J \ ^ SfSj+i -I- satSi I , (1) 

with J < the magnetic coupling strength, N the num- 
ber of magnetic sites, and the spin operator of the jth 
ion with spin s. Their magnetism and spin dynamics 



2 



have intensely been studied by a host of different experi- 
mental techniques'^— the recent experimental work on 
a Crg spin-3/2 ring has probably most clearly demon- 
strated the nature of the elementary excitations.'^'^ 

On the theory side, the rotational-band model (RBM) 
has been established)^""— and there has been striking ex- 
perimental evidence that this picture is a good descrip- 
tion for a variety of finite, bipartite AFM spin clusters 
with varying topologies^""— as well as some frustrated 
systems)^""— but of course also exceptions exist. 

The RBM classifies the lowest-lying magnetic excita- 
tions into two bands with a parabolic energy dependence 
on S, the L and E band {S is the total spin quantum 
number). Physically, the L band can be interpreted as 
rotations of the two AFM sublattices and is, so-to-say, 
a precursor of a long-range Neel-ordered ground state in 
the infinite lattice (it is also known as the tower of states 
or quasi-degenerate joint states^i^). In contrast, the E- 
band states can be associated with quantized spin-wave 
excitations42ii^ For both bands, the energies are propor- 
tional to S{S + 1), but the E band consists of several 
sub-bands which are shifted up in energy by constant 
offsets. For 5 > 1, this can be summarized as 

E{S, k) = \e{^)S{S + 1) + eik) - e(7r), (2) 

with fc = or TT for the L band, where e(fc) may be re- 
garded as the finite-size version of a spin- wave dispersion 
relation, and k would become a wave- vector in an infinite 
system. 

There are several reasons for classifying the low-lying 
excitations into L and E band: the spectrum of a con- 
ventional SWT exhibits zero-energy (Goldstone) modes, 
with wave vector fc = 0, tt, refiecting the (assumed or fac- 
tual) long-range order, which in the context of extended 
magnets are included in the notation of spin waves or 
magnons. However, in a finite system these fc = 0, tt 
modes exhibit a gap with respect to the ground-state en- 
ergy, forming the L band. These excitations are hence 
of different physical significance than the higher-lying 
magnon states [E band). The different nature of the 
L and E band is further manifested by the fact, that 
only the fc = 0,7r modes are responsible for the diverg- 
ing sublattice magnetizations in the conventional SWTs, 
that the L- and £'-band states have different sublattice 
lengths, and that the scaling with N is generally dif- 
ferent. In addition, this distinction is supported by a 
selection rule for the inelastic neutron scattering (INS) 
intensities, which at low temperatures allows only INS 
transitions between L band states or from L band to E 
band states."'^ This key characteristic has been experi- 
mentally demonstrated in great detail)^ Owing to the 
obviously different physical properties in small systems, 
the L band states should be distinguished from the spin 
waves or magnons. 

In the above description, the lowest excited state of an 
AFM Heisenberg ring is the S — \ state of the L band, 
which to a good approximation has a gap of Ac = 4| J|/A^ 




FIG. 1: (Color online) Crystal structure of CsFes [with Fe as 
dark gray spheres (red online) and H atoms omitted]. 

to the S = ground state. Obviously, Ac vanishes for 
N ^ oo; it is thus a property of finite rings. In this con- 
text, Haldane's conjectured^ should also be mentioned. 
It predicts that integer-spin ID chains exhibit a gap in 
the energy spectrum while half-integer ID spin chains do 
not. In small rings, however, there is numerical evidence 
that the gap is largely dominated by Ac without discrim- 
inating integer or half-integer spin4Si^i^ This criterion, 
in fact, may be regarded as a definition of what we mean 
by "small" or "OD" . A further consequence of the small 
size is obviously that the spin waves become quantized, as 
there can be at most as many spin-wave modes as there 
are spin sites in the cluster (the number of spin waves is 
slightly smaller than N, for even AFM rings it is iV — 2). 

In this work we investigate the molecule 
[CsFe8{N(CH2CH20)3}8]Cl, or CsFcg in short (Fig.©. 
By a variety of experimental techniques, such as high- 
field torque magnetometry and INS, the CsFeg molecule 
has been demonstrated to be a good model system for 
an AFM Heisenberg ringi^i^ It is characterized by eight 
Fe(III) ions (s = 5/2), which form an almost planar 
ring with a crystallographic C4 symmetry axis {z axis) 
perpendicular to the wheel plane. The CsFeg molecule 
was shown to be very well described by a generic 
spin Hamiltonian with dominant nearest-neighbor 
AFM Heisenberg couplings plus an uniaxial single-ion 
anisotropy term, 

CN~1 \ N 

^ SjSi+i + satSi ] §1^ . (3) 

i=l / i=l 

All the previous experimental data can consistently be 
described, with an accuracy margin of few percents, with 
the values J = —1.79 meV and D = —0.05 nieV. How- 
ever, all these experiments were only sensitive to the L 
band states. Here we extend the studies to higher ener- 
gies, allowing for an investigation of the ii^-band states. 
The presented INS data will experimentally confirm the 
expected quantized spin-wave modes, complete the pic- 
ture of the excitation modes in the CsFeg molecule, and 
thus establish the only second example of a full experi- 



3 




Energy Transfer (meV) Energy Transfer (meV) 



FIG. 2: INS spectra taken on FOCUS at a temperature of 
T = 1.4 K and wavelengths A = 2.3 A ( main panel) and 3.8 A 
(inset). 



FIG. 3: (Color online) Temperature dependence of the INS 
spectra as recorded on MARI for incident neutron energies 
of 25 meV (main panel) and 15 meV (inset). Intensity was 
summed over all detector banks. 



mental test of the RBM in a small spin cluster. Further- 
more, we will describe the i?-band excitations by different 
SWTs, comparing their performance to each other. This 
provides insight into their merits and drawbacks when 
applied to mesoscopic spin clusters. 

The paper is organized as follows: in the experimental 
Sec. ITT] the INS spectra and S{Q,uj) plots are presented. 
Further, in Sec. IIIII the INS experiments are described 
in terms of the microscopic spin Hamiltonian Eq. ([3]) . In 
Sec. IIVI the i?-band excitations are discussed in a more 
general fashion, revealing their spin-wave character and 
testing different SWTs. Finally, Sec. |V] contains the 
summary and concludes the paper. 



II. EXPERIMENTS 

A powder sample of CsFcs was prepared according 
to literature procedures using non-deuterated starting 
materialsi^ INS data were recorded on the direct time- 
of-flight spectrometer FOCUS at the Paul-Scherrer Insti- 
tute (Villigen, Switzerland) using neutrons with incident 
wavelengths A = 2.3 and 3.8 A, and the MARI spectrom- 
eter at the pulsed neutron source ISIS (Didcot, U.K.) 
using neutrons with incident energies oi Ei = 15 and 
25 meV. The experimental resolutions at the elastic line 
were 1.13 and 0.28 meV for the FOCUS, and 0.39 and 
0.72 meV for the MARI experiments. If not stated oth- 
erwise, the data were summed over all detector banks, 
and the neutron-energy loss side is shown. 

Figure [2] presents the 2.3 A FOCUS data taken at a 
temperature of T = 1.4 K. Clearly, two cold peaks at 
energy transfers of approximately 7 meV (peak III) and 
10 meV (peak IV) are visible. The 3.8 A spectrum, shown 
in the inset of Fig. 1, reveals two additional low-energy 
peaks at energy transfers of about 0.5 meV (peak la) and 
1.3 meV (peak lb). Figure |3] shows INS spectra recorded 
on MARI at temperatures of 5-58 K. The peaks lb. III, 



and IV can again be observed, but thanks to the larger 
energy range and higher resolution, additional features 
become visible. In particular, a further peak at about 
2.3 meV (peak II) is resolved in the 15 meV data (inset 
of Fig.|3]). An inspection of the temperature dependence 
of these four peaks suggests that they are cold transi- 
tions and of magnetic origin. Additional features can be 
observed, such as the broad feature at about 2.5 meV 
(feature marked by "Ph" ) in the high-temperature data 
and the many features at energy transfers above 13 meV. 
Their temperature dependence suggests a phononic ori- 
gin. 

The observed low-energy peaks la, lb, and II are 
fully consistent with the previous investigations on this 
molecule, in particular, the INS experiments-?^'^^ Peaks 
la and lb correspond to the transitions from the S = 
ground state to the first-excited 5* = 1 level, which is 
split by magnetic anisotropy into its components M — 
and Af = ±1 (M is the magnetic quantum number of the 
total spin S). Peak II corresponds to the transition from 
the lowest excited state to the first excited 5* = 2 level. 
The transitions III and IV were hitherto not observed; 
their observation and identification as spin waves form 
the main body of this work. In the following, we will 
analyze the data more carefully in order to demonstrate 
the magnetic nature of the peaks I to IV, and the lattice 
origin of the other features in the spectra. In doing so, we 
will focus on the two peaks III and IV (as the magnetic 
peaks at lower energies are already well understood). 

First, we analyze the full S{Q,uj) dependence as pre- 
sented in Fig. m for the temperatures 5 and 58 K and 
settings Ei = 15 and 25 meV. The Q dependence of 
the scattering intensity allows a clear differentiation: the 
magnetic scattering is determined by the magnetic form 
factors of the involved metal ions and interference factors 
refiecting their spatial arrangement,— it is hence most 
significant at lower Q values. In contrast, the phonon 
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FIG. 4: (Color online) S{Q,w) plots for the temperatures T — 5 K (upper plots) and 58 K (lower plots) for the two incident 
neutron energies Ei — 15 meV (left) and 25 meV (right) as measured on MARL The intensity is color-coded (blue = low, red 
— high intensity); the color scale is identical for each of the four plots. 



scattering intensity increases with Q as Q^. In the low- 
temperature data, peaks III and IV (marked by black ar- 
rows) are visible in the whole Q range, in particular, at 
small Q values, which strongly supports a magnetic origin 
of these peaks. In the high Q regime the phonon-induced 
scattering clearly dominates. The phonon scattering also 
increases strongly with temperature, as demonstrated by 
the 58 K data, where scattering is dominated by phonons 
at all momentum transfers. 

Figure [5] shows cuts along the energy-axis of the 5 K 
spectra for three different regimes of momentum transfer 
(black, red, and green circles). It reveals that the inten- 
sity of the peaks lb, II, III, and IV roughly stays constant 
with Q. It might be surprising that these peaks exhibit 
significant intensity even at Q values as large as approxi- 
mately 4 A, where the magnetic form factors have already 
dropped considerably; we attribute this to a strong inco- 
herent scattering in our non-deuterated sample smearing 
out the Q dependence. These plots further demonstrate 
that the feature at approximately 2.5 meV (compare with 
Fig. [2|) must indeed be a phonon. 

In Fig. [Sj the gray curves were obtained by scaling the 
T = 58 K data by the Bose factor [l-exp{-hxj/kBT)]-\ 
which governs the temperature dependence of phonon 
scatteringj^ Considering that at high temperatures the 
magnetic scattering intensity is strongly reduced as it 
spreads out over essentially all energies, this curve should 
provide an estimate of the phonon background. A com- 
parison with the 5 K data nicely confirms the magnetic 
nature of peaks III and IV, and also proves that the 
scattering at higher energies above 12 meV is purely 



phononic. The data after correction by the Bose factor 
is shown in Fig. [SJa) as a S{Q,uj) plot for Ei = 25 meV, 
T = 5 K. The two peaks III and IV are clearly visible, 
proving their magnetic origin. 

The linewidths of the low-energy peaks la, lb, and II lie 
within experimental accuracy determined by the instru- 
mental resolution. In contrast, the linewidths of peaks 
III and IV are considerably broadened, as is most clearly 
seen in Fig. ^h). However, our analysis will reveal that 
they consist of several close-lying transitions. Finally, 
it is interesting to note that in the energy regime above 
peak IV, i.e., above approximately 13 meV, no significant 
magnetic scattering intensity could be observed. In this 
regime, the scattering is entirely of phononic origin, as is 
evident from Fig. |6l^a) and even more so from Fig.[5][a). 



III. ANALYSIS 

We have simulated the INS spectrum of Eq. ([3]) by 
calculating the energies and wave functions via a sparse- 
matrix exact numerical diagonalization method and us- 
ing the formulas of Refs. 68 and J.O. In order to find 
the parameters which yield the best fit to the exper- 
imental spectra we have systematically scanned the J 
and D parameter space, yielding J = —1.79(5) meV and 
D = —0.05(7) meV. These values are in excellent agree- 
ment with the previous findings.— The dependence of 
the simulated INS spectrum on J and D is shown in 
Figs.EJa) andlT^b), respectively. The low-energy part of 
the spectrum is affected by both J and D. In contrast. 
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FIG. 5: (Color online) Q slices for (a) the Ei = 25 meV and 
(b) the Ei = 15 meV data at T = 5 K. In panel (a), the lowest 
(black), middle (red), and upper (green) curves correspond to 
Q e [1.0, 2. 5[, [2.5,4.0[, and [4.0,5.5] A-\ and in panel (b) to 
Q € [1.0, 2. 0[, [2.0,3.0[, and [3.0,4.0] A'^ The curves have 
been offset in order to improve visibility. The gray curves 
represent the Bose-factor corrected 58 K data as discussed in 
the text. The solid lines at the bottom are calculated INS 
spectra obtained with the best-fit parameters. The spikes 
(red) mark the exact position and relative scattering strengths 
of the individual INS transitions, the solid (black) curve was 
obtained by convolution with a Gaussian accounting for the 
experimental resolution. 




Momentum Transfer (A ^) 

FIG. 6: (Color online) (a) S{Q,lo) plot of the Ei = 25 meV 
MARI data at T = 5 K after subtraction of the phonon back- 
ground as described in the text, (b) S{Q, uj) plot (T = 5 K) as 
calculated from the spin Hamiltonian Q with the parameters 
J = —1.79 meV and D — —0.05 meV (a Gaussian linewidth 
of 0.72 meV was used). 



the high-energy part is predominantly affected by J. In- 
deed, a 10% variation in J has a much larger effect than 
a 50% variation in D. This explains the rather large es- 
timated standard deviation of the D value that we have 
determined. The simulated INS spectra corresponding to 
the best-fit values are shown at the bottom of Figs. EJa) 
and[5]Jb) as solid lines, and as a S{Q, uj) plot in Fig.lHIb); 
the good agreement with the data is obvious. 

Figure [5] shows the relevant part of the calculated spec- 
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FIG. 7: (Color online) (a) Dependence of the INS spectrum 
on J, which is varied by a maximum of ±10% around a center 
value of J = —1.79 meV (red curve). The anisotropy is D = 
—0.05 meV. (b) Dependence of the INS spectrum on D, which 
is varied by a maximum ±50% around a center value of D = 
—0.05 meV (red curve). The coupling is J = —1.79 meV. (c) 
Dependence of the INS spectrum on a modulation of the J 
values along the ring. For all curves the anisotropy is D = 
—0.05 meV. Red (solid) curve: no modulation, J\ = J2 = 
— 1.79 meV. Blue (short dashed) curve: Ji = —2.14 meV and 
J2 = —1.45 meV. Black (long dashed) curve: Ji = —2.48 meV 
and J2 = —1.10 meV. 




Total Spin S 

FIG. 8: (Color online) Calculated low-energy spectrum of 
Hamiltonian Q with J = -1.79 meV and D = -0.05 meV. 
Only states in the sector S <2 are plotted. The energy of the 
ground state was set to zero. Due to the anisotropy term, the 
degeneracy in M of each S multiplet is partially lifted. The 
relevant INS transitions are plotted as solid green (dashed 
red) arrows, marking transitions to the E band (L band). 
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TABLE I: Calculated contributions to the observed INS tran- 
sitions. The quantum states are labeled as \S,M,k), as in 
Fig. El The ground state is lgs> = |0,0,0>. \i) and |/) de- 
notes the initial and final state in the transition. Ef — Ei is 
the excitation energy in units of meV, and hf indicates the 
transition strength, normalized to the strongest transition la. 
Transitions with strength hf < 0.05 have been omitted. 







f \ 






la 


|gs) 


|l,0,7r) 


0.50 


1.00 


lb 


|gs) 


|1,±1,^> 


1.34 


0.97 


II 


|l,±l,7r) 


12, ±2,0) 


2.05 


0.06 




|l,0,7r) 


|2,0,0) 


2.25 


0.12 




|1,0,^) 


|2,±1,0) 


2.30 


0.20 


III 


|l,0,7r) 


|l,0,37r/4) 


6.83 


0.09 




|gs) 


|l,±l,37r/4) 


7.30 


0.33 




|gs) 


|1,±1,V4> 


7.33 


0.05 




Igs) 


|l,0,37r/4) 


7.61 


0.13 


IV 


|gs) 


|1,0,V2> 


10.59 


0.05 




|gs) 


|1,±1,V2> 


10.15 


0.14 



trum, detailing the assignment of the observed peaks. 
The assignment of the lower- ene rgy peaks has already 
been discussed in Refs. IH and HI and in the above. 
The peaks III and IV correspond to transitions from the 
S = ground state to higher-lying 5 = 1 levels, which 
are slightly split due to the magnetic anisotropy. Accord- 
ingly, each of the peaks consists of a number of close- lying 
transitions, giving rise to the experimentally observed 
broadening of the peaks. Furthermore, at the measure- 
ment temperature of 5 K, also the transitions from the 
first excited state weakly contribute. The INS intensities 
of the individual transitions are indicated by the verti- 
cal thin lines in Fig. [5] The calculated contributions to 
all observed INS peaks, i.e., the quantum numbers of 
the involved states, the transition energies, and the scat- 
tering intensities, are detailed in Table ID The model 
Hamiltonian ([3]) assumes a perfect symmetry of the ring 
(i.e., an eight-fold symmetry axis). However, though it 
is close to Ss, the molecular symmetry of CsFeg is C4. 
Hence, weak deviations exist raising the question of how 
important they are for the magnetism in CsFcs- This 
point is relevant for the understanding and explanation 
of phenomena such as the field-induced spin- Jahn- Teller 
effecli^"— or the significance of possible Dzyaloshinski- 
Moriya (DM) interactions)^ In this context, the present 
experiment is not sensitive to such small effects as a DM 
interaction or a biaxial anisotropy term {E term), but 
a variation in the exchange-coupling constants along the 
ring, e.g., in a J1-J2 alternating fashion, can be tested. 
Indeed, from general arguments it is expected that the 
L band states are rather insensitive to such a variation 
(as their energies are affected only in second order, see 
Ref . l65h . However, this is not the case for the E band or 
spin-wave states, respectively, which hence provide a bet- 
ter means for a check. A simulation of the INS spectrum 



at 5 K with a modulation of the exchange-coupling con- 
stant of up to ±40% around its best-fit value is shown 
in Fig. As expected, the low-energy part is not 

much affected but the spin-wave excitations at 7 and 
10 mcV have significantly moved toward lower energies. 
By comparison with the experimental widths of the peaks 
we conclude that a modulation of the exchange-coupling 
constants is smaller than ±20% in CsFeg. 



IV. DISCUSSION 

In the previous sections the following points have been 
made: peaks HI and IV are of magnetic origin, they are 
well reproduced by the generic spin Hamiltonian ([3]) , they 
correspond to transitions from the ground state to ex- 
cited S = 1 states, and there is negligible magnetic scat- 
tering at energies above peak IV. In this section we will 
discuss some implications of these findings. 

As presented in Sec. [H the low-energy spectrum of 
AFM Heisenberg spin rings can be described by the 
REM, which subdivides the spectrum into the L and 
E band; all other states are denoted collectively as the 
quasi-continuum. The L band is formed by the lowest 
spin multiplet in each sector of total spin S. For CsFes, 
the transitions within the states of the L band were al- 
ready investigated intensely beforej^iS in the present 
work they are observed as peaks la, lb, and II. The E 
band consists of the N — 2 next higher-lying spin multi- 
plets in each spin sector S* > 1. These states are exactly 
those which are observed as the two peaks HI and IV 
in our INS data. Only two broad INS peaks appear be- 
cause of the (near) degeneracy of some states, which are 
due to the ring symmetry and the two-sublattice struc- 
ture in the AFM ring, as will also become clear in the 
following. The REM also predicts a selection rule for 
the matrix elements, which in fact is a key aspect of the 
REM as it ultimately justifies the distinction of the states 
into the L band, E band, and quasi-continuum. It states 
that, starting from the ground state, only transitions to 
states of the L or E band have significant INS intensities, 
while transitions into the quasi-continuum are negligibly 
weak. Thus, at low temperatures no magnetic INS in- 
tensity should be observed at energies above the E band 
or peak IV, respectively, exactly as observed in our ex- 
periment. Therefore, in conjunction with the previous 
work j^^i^'' which confirmed the L band in great detail, 
the data presented here provides a thorough confirmation 
of the REM for CsFeg. This is only the second example 
where this has been achieved (the first being on the Crg 
molecule, Ref. [soh . 

The connection between E band and quantized AFM 
spin waves has been discussed before in the literature to 
some extent In the following it will be made 
explicit by applying various SWTs available in the lit- 
erature to a small AFM Heisenberg ring. The analysis 
will also give insight into the strengths and weaknesses 
of these methods when applied to small spin clusters. 
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As it was shown in Sec. Ill, the energies of peaks III 
and IV do not depend much on the magnetic anisotropy. 
Hence, it is ignored in the calculations in this section, i.e., 
Eq. ([T]) is used. Excitation energies are given relative to 
the ground-state energy Eg. 

It is useful to first discuss the symmetries of the model 
and Fig. [8l The AFM Heisenberg ring exhibits spin- 
rotational symmetry, which gives rise to quantum num- 
bers S and M, and a. Dn spatial symmetry. In a pure 
Heisenberg model the states are degenerate with respect 
to M. In CsFcg this degeneracy is lifted by the mag- 
netic anisotropy, as indicated in Fig. |8l but the splitting 
is weak for the higher-lying states, and a classification 
in terms of spin multiplets or S", respectively, is useful. 
As regards the spatial symmetry, it is convenient to re- 
lax it to Cat. The quantum numbers, which we will call 
shift quantum number are related to a shift operator 
f\q) = exp(ig27r/iV)|g) with q = 0,l, ...,N ~ 1. Alterna- 
tively, one may use "wave vectors" k = 2nq/N + const, 
where we chose const such that for the ground state 
k = Q. The actual symmetry, however, is D^, which 
requires that states with q and N — q, or k and — fc, re- 
spectively, are degenerate. We now turn to Fig. [H The 
S' = sector embraces the ground state with k = 0. The 
remaining S = levels are high up in energy and are 
part of the quasi-continuum (not shown in Fig. In 
the 5* = 1 sector, the (three) lowest states correspond to 
the first-excited spin multiplet of the L band, for which 
k = n. This multiplet appears split in Fig. [S] because of 
the magnetic anisotropy. The next iV — 2 = 6 higher- 
lying spin multiplets with k = ±^,±-|,±^ form the E 
band. Their counting is difficult in Fig. [8] because of the 
anisotropy splitting and the degeneracies. This pattern 
of states continues for S* > 1 (not shown in Fig. |8]). 

Alternatively to the energy-vs-S* representation, the 
structure of the energy spectrum may be discussed in an 
energy-vs-fc representation, which for an infinite system 
would be the dispersion relation. However, in a small 
system, with no periodicity along an extended axis, k 
in general does not correspond to a "real" wave vector 
(although it can approach one for iV — ^ c» in the case 
of a ring). Furthermore, k may assume only a set of fi- 
nite values, and in consequence, the allowed energies are 
discrete values, yielding quantized spin waves in our con- 
text. All available SWTs can easily be applied to small 
spin clusters if one chooses to work in real space and not 
momentum space (see Refs.lsolandlTSl). Generally, SWTs 
work with bosonic excitations from the ground state, the 
magnons or spin waves. The one-magnon excitations cor- 
respond to the L- and i?-band states in the S = 1 sector. 
Other states may, in principle, be generated by multiple- 
magnon excitations, but because of the technical diffi- 
culties SWT is rarely developed to this stage. In other 
words, the SWTs provide us with the energies of the L- 
and i5-band states in the S = 1 sector or the values e{k) 
in Eq. respectively. The dependence on S is essen- 
tially out of their reach. The excitation energies of these 
states are plotted as a function of k in Fig.lHlfor different 



SWTs as well as the exact numerical result. 

In the early 1950s, the (conventional) linear and in- 
teracting SWTs were establishedJ^ii^ Low-dimensional 
systems present a challenge, since divergencies appear in 
ID (and OD) due to the fact that in these SWTs the 
modes with k = and tt are degenerate and have zero 
energy (Goldstone modes). The L band is related to 
these modes; hence the gap between the S — ground 
state and the first-excited S = 1 state (singlet-triplet 
gap), which in small AFM systems is always present, 
is zero, i.e., can never be reproduced by these SWTs. 
Later on, the modified SWTsii"— and Schwinger-boson 
mean-field theory (SBMFT)^i22 were introduced, which 
eliminate the drawback of a zero singlet-triplet gap, sug- 
gesting their applicability to small spin clusters. Fortu- 
nately, for all these theories, linear SWT (LSWT), inter- 
acting SWT (ISWT), hnear modified SWT (LMSWT), 
interacting modified SWT (IMSWT), and SBMFT, ana- 
lytical results are available for ID, which can be directly 
applied to the AFM ring (various IMSWTs exist, we have 
used the fuU-diagonalization IMSWT of Ref. [Zi). Inter- 
estingly, the SBMFT yields exactly the same excitation 
energies as the LMSWT, and also the ground-state en- 
ergy agree after introduction of a correction factor i ^^i^^ 
Hence, there is no need to consider the SBMFT sepa- 
rately in the following. The singlet-triplet gap may be 
approximated by Ac = 4| Jj/A'^, which becomes exact in 
the classical limit J^i^i^ Adding this gap to the exci- 
tation energies of the conventional SWTs is tempting 
(and may be justified, see Refs. i21| and i23i) . but does 
not improve the ISWT; the result for LSWT, denoted as 
LSWT -I- Ac, is discussed below. 

The excitation energies as function of k are presented 
in Fig. [9l The interpretation of the INS transitions in 
terms of spin waves is obvious: peak IV at 10 meV corre- 
sponds to the one-magnon excitations with /c = and 
peak HI at 7 meV is related to excitations of magnons 
with k = ±j and The "quantization" of the spin- 

wave energies due to the system's finite size can also be 
regarded as being due to a "confinement" of the magnons 
in three-dimensional space. The dispersion relation as 
calculated by SWT is also symmetric with respect to 
k = ±j due to the AFM sublattice structure. 

Table II compiles the ground-state energy Eg, singlet- 
triplet gap A,— width e(7r/2) — e(7r) of the excita- 
tion spectrum, and mean deviation = Sfe[^(^) ~ 
£c^act{k)]'^ / {^N\J\s) from the exact energies as calcu- 
lated by exact numerical diagonalization and the various 
SWTs for the case iV = 8, s = 5/2. Inspection of Table II 
and Fig. [9] shows: 

(1) As discussed beforehand the conventional SWTs 
(LSWT and ISWT) do not produce a singlet-triplet gap 
while the modified SWTs do. Yet, they grossly underes- 
timate the gap, by a factor of 2. In contrast, the clas- 
sical estimate Ac is off by only 7%. Hence, although 
the modified SWTs reproduce the singlet-triplet gap, 
they do not provide a good description of the low-energy 
or low-temperature behavior, respectively. Accordingly, 
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FIG. 9: (Color online) Spin-wave excitation energies as cal- 
culated by exact numerical diagonalization (stars) normalized 
to I J| and the different SWTs (solid, dashed and dotted lines) 
as discussed in the text. The abbreviations are LSWT: linear 
SWT, ISWT: interacting SWT, IMSWT: fuU-diagonalization 
interacting modified SWT, LSWT + A^: linear SWT plus 
classical gap Ac. 

TABLE II: Exact and SWT results for the singlet-triplet gap, 
energy width, and mean deviation for a A'^ = 8, s = 5/2 ring. 

EJ\J\ A/|J| Kfc = 7r/2) - A]/| J| 



Exact -58.1105 0.537 5.180 

LSWT -57.2600 5 0.133 

LMSWT -57.8049 0.248 4.758 0.115 

ISWT -58.2433 5.396 0.093 

IMSWT -58.1095 0.268 5.129 0.065 

LSWT-f-A, 0.5 5 0.025 



estimate the width. It is noted that the modified SWTs 
generally predict a smaller width than the corresponding 
conventional SWTs. Further, among the SWTs, IMSWT 
yields the best assessment of the excitations. Notably, 
the empirical LSWT-f-Ac calculation performs even bet- 
ter than the IMSWT. 

(3) Regarding the ground-state energy, the IMSWT 
yields the best results, with a remarkably small deviation 
of only 0.002%. This implies that also the ground-state 
wave function is accurately obtained by this method. It 
would be interesting to verify this in the future. 



V. CONCLUSION 

In conclusion, we have observed by INS the excitation 
of quantized spin- wave modes in the molecule CsFcg, i.e., 
an eight-membered AFM Heisenberg ring. Only a small 
number of discrete points on the e{k) spin-wave disper- 
sion curve are allowed due to the confinement of the mag- 
netic excitations to the ring. Further, the INS data re- 
vealed that at low temperatures the magnetic scattering 
intensity at energies above the most energetic spin- wave 
excitation is negligibly small. These findings, combined 
with the demonstration of the L band in previous works, 
establish a thorough confirmation of the RBM for CsFeg. 
We have compared several SWTs with numerically ex- 
act results. The fuU-diagonalization interacting modified 
SWT delivers excellent approximations to the ground- 
state energy and spin-wave excitations, but the singlet- 
triplet gap and hence the low-temperature properties of 
the ring are only qualitatively reproduced. Regarding the 
excitation energies the empirical LSWT -I- Ac calculation 
performs best. 



LSWT-I-Ac performs best. 

(2) The width e(7r/2) — e(7r) of the spin- wave excitation 
spectrum is obtained reasonably well in all considered 
SWTs, with the expected trend that the interacting ver- 
sions do better. Except for the ISWT, the SWTs under- 



Acknowledgments 

The authors thankfully acknowledge funding by the 
Deutsche Forschungsgemeinschaft. 



* Electronic address: |jan.drei ser@physik.uni-freiburg.de| 

t Electronic address: |oliver. w aldmann@physik.uni-freiburg.de 

* Present address: Spallation Neutron Source, Oak Ridge 
National Laboratory, Oak Ridge, Tennessee 37831, USA 

^ A. R. Rocha, V. M. Garcia-suarez, S. W. Bailey, C. J. 
Lambert, J. Ferrer, and S. Sanvito, Nature Mater. 4, 335 
(2005). 

^ L. Bogani and W. Wernsdorfer, Nature Mater. 7, 179 
(2008). 

^ C. Schlegel, J. van Slageren, M. Manoli, E. K. Brechin, 
and M. Dressel, Phys. Rev. Lett. 101, 147203 (2008). 

* S. Bertaina, S. Gambarelli, T. Mitra, B. Tsukerblat, 
A. Miiller, and B. Barbara, Nature 453, 203 (2008). 

^ M. Mannini, F. Pineider, P. Sainctavit, C. Danieli, 
E. Otero, C. Sciancalepore, A. M. Talarico, M.-A. Arrio, 



A. Cornia, D. Gatteschi, and R. Sessoli, Nature Mater. 8, 
194 (2009). 

® J. des Cloizeaux and J. J. Pearson, Phys. Rev. 128, 2131 
(1962). 

^ Zheng Weihong and C. J. Hamer, Phys. Rev. B 47, 7961 
(1993). 

* P. V. Hendriksen, S. Linderoth, and P.-A. Lindgard, Phys. 

Rev. B 48, 7259 (1993). 
^ C. Mathieu, J. Jorzick, A. Frank, S. O. Demokritov, A. N. 

Slavin, B. Hillebrands, B. Bartenlian, C. Chappert, D. De- 

canini, F. Rousseaux, and E. Cambril, Phys. Rev. Lett. 81, 

3968 (1998). 

" J. Jorzick, S. O. Demokritov, B. Hillebrands, M. Bailleul, 
C. Fermon, K. Y. Guslienko, A. N. Slavin, D. V. Berkov, 
and N. L. Gorn, Phys. Rev. Lett. 88, 047204 (2002). 



9 



J. Podbiclski, F. Giesen, and D. Grundler, Phys. Rev. Lett. 

96, 167207 (2006). 

J. Topp, J. Podbiclski, D. Heitmann, and D. Grundler, 

Phys. Rev. B 78, 024431 (2008). 
^■^ P. W. Anderson, Phys. Rev. 86, 694 (1952). 
" R. Kubo, Phys. Rev. 87, 568 (1952). 

T. Oguchi, Phys. Rev. 117, 117 (1960). 
^® M. Takahashi, Phys. Rev. Lett. 58, 168 (1987). 
" M. Takahashi, Phys. Rev. B 40, 2494 (1989). 

,J. E. Hirsch and S. Tang, Phys. Rev. B 40, 4769 (1989). 

S. Tang, M. E. Lazzouni, and J. E. Hirsch, Phys. Rev. B 

40, 5000 (1989). 
^° S. M. Rczcndc, Phys. Rev. B 42, 2589 (1990). 
^1 Q. F. Zhong and S. Sorella, Europhys. Lett. 21, 629 (1993). 

C. Lavalle, S. Sorella, and A. Parola, Phys. Rev. Lett. 80, 

1746 (1998). 

A. E. Trumper, L. Capriotti, and S. Sorella, Phys. Rev. B 
61, 11529 (2000). 

S. Itch, Y. Endoh, K. Kakurai, and H. Tanaka, Phys. Rev. 
Lett. 74, 2375 (1995). 

N. E. Ivanov and P. C. Ivanov, Phys. Rev. B 46, 8206 
(1992). 

^® R. Wieser, E. Y. Vedmedenko, and R. Wiesendanger, Phys. 

Rev. Lett. 101, 177202 (2008). 
'^^ S. T. Ochsenbein, O. Waldmann, A. Sieber, G. Carver, 

R. Bircher, H. U. Giidel, R. S. G. Davies, G. A. Timco, 

R. E. P. Winpenny, H. Mutka, and F. Fernandez-Alonso, 

Europhys. Lett. 79, 17003 (2007). 
^® S. T. Ochsenbein, F. Tuna, M. Rancan, R. S. G. Davies, 

C. A. Muryn, O. Waldmann, R. Bircher, A. Sieber, 
G. Carver, H. Mutka, F. Fernandez-Alonso, A. Podlesnyak, 
L. P. Engelhardt, G. A. Timco, H. U. Giidel, and R. E. P. 
Winpenny, Chem. Eur. J. 14, 5144 (2008). 

G. Chaboussant, A. Sieber, S. Ochsenbein, H.-U. Giidel, 
M. Murrie, A. Honecker, N. Fukushima, and B. Normand, 
Phys. Rev. B 70, 104422 (2004). 

^° O. Waldmann, T. Guidi, S. Carretta, C. Mondelli, and 

A. L. Dearden, Phys. Rev. Lett. 91, 237202 (2003). 

S. Yamamoto and T. Nakanishi, Phys. Rev. Lett. 89, 
157603 (2002). 

H. Hori and S. Yamamoto, Phys. Rev. B 68, 054409 (2003). 

D. Gatteschi and R. Sessoli, Molecular Nanomagnets (Ox- 
ford University Press, 2006). 

B. Barbara, L. Thomas, F. Lionti, I. Chiorescu, and 
A. Sulpice, J. Magn. Magn. Mater. 200, 167 (1999). 

G. D. Christou, and D. N. Hendrickson, MRS Bull. 25, 66 
(2000). 

^® D. Gatteschi and R. SessoU, Angew. Chem. Int. Ed. 42, 
268 (2003). 

W. Wernsdorfer, Adv. Chem. Phys. 118, 99 (2007). 
^® K. L. Taft, C. D. Delfs, G. C. Papaefthymiou, S. Foner, 
D. Gatteschi, and S. J. Lippard, J. Am. Chem. Soc. 116, 
823 (1994). 

M.-H. Julien, Z. H. Jang, A. Lascialfari, F. Borsa, M. Hor- 
vatic, A. Caneschi, and D. Gatteschi, Phys. Rev. Lett. 83, 
227 (1999). 

A. Cornia, A. G. M. Jansen, and M. Affronte, Phys. Rev. 
B 60, 12177 (1999). 

,J. van Slageren, R. Sessoh, D. Gatteschi, A. A. Smith, 
M. HeUiwell, R. E. P. Winpenny, A. Cornia, A. L. Barra, 
A. G. M. Jansen, E. Rentschler, and G. A. Timco, Chem. 
Eur. J. 8, 277 (2002). 

S. Carretta, J. van Slageren, T. Guidi, E. Liviotti, C. Mon- 
delli, D. Rovai, A. Cornia, A. L. Dearden, F. Carsughi, 



M. Alfronte, C. D. Frost, R. E. P. Winpenny, D. Gat- 
teschi, G. Amoretti, and R. Caciuffo, Phys. Rev. B 67, 
094405 (2003). 

B. Pilawa, R. Boffinger, I. Keilhauer, R. Leppin, I. Oden- 
wald, W. Wendl, C. Berthier, and M. Horvatic, Phys. Rev. 
B 71, 184419 (2005). 

P. W. Anderson, Basic Notions of Condensed Matter 
Physics (Bonjamin/Cummings, Menlo Park, 1984). 
*^ B. Bernu, C. Llmillier, and L. Pierre, Phys. Rev. Lett. 69 
2590 (1992). 

*® B. Bernu, P. Lecheminant, C. Lhuillier, and L. Pierre 
Phys. Rev. B 50, 10048 (1994). 

A. Chiolero and D. Loss, Phys. Rev. Lett. 80, 169 (1998) 
*® J. Schnack and M. Luban, Phys. Rev. B 63, 014418 (2000) 

O. Waldmann, Phys. Rev. B 65, 024424 (2001). 
^° C. Lhuillier, arXiv:cond-mat/0502464 (2005). 

T. Guidi, S. Carretta, P. Santini, E. Liviotti, N. Magnani 

C. Mondelli, O. Waldmann, L. K. Thompson, L. Zhao 

C. D. Frost, G. Amoretti, and R. Caciuffo, Phys. Rev. B 
69, 104432 (2004). 

R. Caciuffo, T. Guidi, G. Amoretti, S. Carretta, E. Liv- 
iotti, P. Santini, C. MondeUi, G. Timco, C. A. Muryn, and 
R. E. P. Winpenny, Phys. Rev. B 71, 174407 (2005). 
O. Waldmann, Coord. Chem. Rev. 249, 2550 (2005). 
P. Santini, S. Carretta, G. Amoretti, T. Guidi, R. Caciuffo, 
A. Caneschi, D. Rovai, Y. Qiu, and J. R. D. Copley, Phys. 
Rev. B 71, 184405 (2005). 

O. Waldmann, H. U. Giidel, T. L. Kelly, and L. K. Thomp- 
son, Inorg. Chem. 45, 3295 (2006). 
^® T. Guidi, J. R. D. Copley, Y. Qiu, S. Carretta, P. Santini, 
G. Amoretti, G. Timco, R. E. P. Winpenny, C. L. Dennis, 
and R. Caciuffo, Phys. Rev. B 75, 014408 (2007). 
J. Schnack, M. Luban, and R. Modler, Europhys. Lett. 56, 
863 (2001). 

^® V. O. Garlea, S. E. Nagler, J. L. Zarestky, C. Stassis, 

D. Vaknin, P. Kogerler, D. F. McMorrow, C. Niedermayer, 

D. A. Tennant, B. Lake, Y. Qiu, M. Exler, J. Schnack, and 
M. Luban, Phys. Rev. B 73, 024414 (2006). 

^® O. Waldmann, Phys. Rev. B 75, 012415 (2007). 

^° G. Miiller, H. Thomas, H. Bock, and J. C. Bonner, Phys. 

Rev. B 24, 1429 (1981). 
®^ F. D. M. Haldane, Phys. Rev. Lett. 50, 1153 (1983). 
®2 J. Schnack, Phys. Rev. B 62, 14855 (2000). 
®^ L. Engelhardt and M. Luban, Phys. Rev. B 73, 054430 

(2006). 

®* R. W. Saalfrank, I. Bernt, E. UUer, and F. Hampel, Angew. 

Chem., Int. Ed. Engl. 36, 2482 (1997). 
®® O. Waldmann, R. Koch, S. Schromm, J. Schiilein, 

P. Miiller, I. Bernt, R. Saalfrank. F. Hampel, and 

E. Balthes, Inorg. Chem. 40, 2986 (2001). 

O. Waldmann, C. Dobe, H. U. Giidel, and H. Mutka, Phys. 

Rev. B 74, 054429 (2006). 
^'^ O. Waldmann, C. Dobe, H. Mutka, A. Furrer, and H. U. 

Giidel, Phys. Rev. Lett. 95, 057202 (2005). 
®* O. Waldmann, Phys. Rev. B 68, 174406 (2003). 
®^ G. Shirane, S. M. Shapiro, and J. M. Tranquada, Neutron 

Scattering with a Triple-Axis Spectrometer (Cambridge 

University Press, Cambridge, England, 2002). 
™ O. Waldmann and H. U. Giidel, Phys. Rev. B 72, 094422 

(2005). 

O. Waldmann, C. Dobe, S. T. Ochsenbein, H. U. Giidel, 
and I. Sheikin, Phys. Rev. Lett. 96, 027206 (2006). 
''^ A. Soncini and L. F. Chibotaru, Phys. Rev. Lett. 99, 
077204 (2007). 



10 



V. Lante, I. Rousochatzakis, K. Pcnc, O. Waldmann, and 
F. Mila, Phys. Rev. B 79, 180412(R) (2009). 
F. Cinti, M. Affronte, and A. Jansen, Eur. Phys. J. B 30, 
461 (2002). 

O. Cepas and T. Ziman, Prog. Theor. Phys. Suppl. 159, 
280 (2005). 

A. Auerbach and D. P. Arovas, Phys. Rev. Lett. 61, 617 
(1988). 

S. Sarker, C. Jayaprakash, H. R. Krishnamurthy, and 



M. Ma, Phys. Rev. B 40, 5028 (1989). 

S. Yamamoto and H. Hori, J. Phys. Soc. Jpn. 72, 769 

(2003). 

D. P. Arovas and A. Auerbach, Phys. Rev. B 38, 316 
(1988). 

We define the singlet-triplet gap in the case of SWTs as 
A = e(7r). 



